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Abstract 

The electric field used in most electrorheological (ER) experiments is usually 
quite high, and nonlinear ER effects have been theoretically predicted and ex- 
perimentally measured recently. A direct method of measuring the nonlinear 
ER effects is to examine the frequency dependence of the same effects. For a 
sinusoidal applied field, we calculate the ac response which generally includes 
higher harmonics. In is work, we develop a multiple image formula, and cal- 
culate the total dipole moments of a pair of dielectric spheres, embedded in a 
nonlinear host. The higher harmonics due to the nonlinearity are calculated 
systematically. 

PACS Number(s): 83.80.Gv, 72.20.Ht, 82.70.Dd, 41.20.-q 
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I. INTRODUCTION 



Electrorheological (ER) fluids consist of liiglily polarizable particles in a nearly insulat- 
ing fluid. Upon the application of electric fields, the apparent viscosity of ER fluids can 
be changed by several orders of magnitude, due to the formation of chains of and columns 
of particles across the electrodes in the direction of the applied field. The rapid transi- 
tion between the fluid and solid phases renders this material the potential of important 
technological applications. 

On the other hand, the applied electric field used in most ER experiments is usually 
quite high, and important data on nonlinear ER effects induced by a strong electric field 
have been reported by Klingenberg and coworkers Recently, the effect of a nonlinear 
characteristic on the interparticle force has been analyzed in an ER suspension of nonlinear 
particles and further extended to a nonlinear host medium 0. These work confirmed 
the previous theoretical results that the attractive force between two touching spheres can 
have a quasi-linear dependence 

A convenient method of probing the nonlinear characteristic is to measure the harmonics 
of the induced polarization under the application of a sinusoidal (ac) electric field . When 
a nonlinear composite with nonlinear dielectric particles embedded in a host medium, or 
with a nonlinear host medium is subjected to a sinusoidal field, the electrical response in the 
composite will in general be a superposition of many sinusoidal functions P] . It is natural to 
investigate the effects of a nonlinear characteristic on the interparticle force in an ER fluid 
which can be regarded as a nonlinear composite medium The strength of the nonlinear 
polarization is reflected in the magnitude of the harmonics. 

In this work, we will develop a self-consistent theory to calculate the ac response of a 
nonlinear ER fluid. Our theory goes beyond the simple point-dipole approximation and 
accounts for the multipole interaction between the polarized particles. 
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II. NONLINEAR POLARIZATION AND ITS HIGHER HARMONICS 



We first consider an isolated spherical particle of radius a and dielectric constant e.p. 
The sphere is placed in a host medium of dielectric constant e^- When a time-dependent 
electric field E = ^'(t)^ is applied. The particle will be polarized and its surface charge will 
contribute to it a dipole moment pq = emo!^ (3E{t), where {3 is the dipolar factor given by: 
/3 = (ep — em)/(ep + 2em)- In ER fluids, when the particles get close to one another, they 
do not behave as point dipole and we have to consider the mutual polarization between the 
particles. 

Next, we examine the effect of a nonlinear characteristic on the induced dipole moment. 
We concentrate on the case where the host medium have a nonlinear dielectric constant, while 
the suspended particles are linear. The nonlinear characteristic gives rise to a field-dependent 
dielectric coefficient. In other words, the electric displacement-electric field relation in the 
host medium is given by: 

where Xm is the nonlinear coefficient of the particles. We have assumed that the average 
field inside the host medium is uniform. This is called the decoupling approximation p|. It 
has been shown that such an approximation yields a lower bound for the accurate result . 
Furthermore, in ER fluids the frequency of the applied field is low; we can assume both e and 
X are independent of frequency. As a result, the induced dipole moment is field-dependent 
and is given by: 

Po = e^a^^Eit), (2) 
where (5 is the field-dependent dipolar factor and is given by: 

^~ ep + 2em~ ep + 2{em + Xm{El))' ^ ' 

When the polarized spheres approach to one another, they will be further polarized by 
the mutual polarization effect. As a result, the point-dipole approximation breaks down and 
we must consider the multipole moments. 



Let us consider a pair of nonlinear dielectric spheres of the same radius a, separated by 
a distance r. Each of them has a dielectric coefficient e^. By using the method of multiple 
images [p|,p!Q[, we deduce the total dipole moment of the spheres: 
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Pr-P.£(-.grf .r''°, I ■ (4) 



n=0 

The subscript T denotes that the applied electric field is perpendicular to the line joining 
the centers of the particles (i.e., a transverse field). The parameter a is related to the 
separation between the particles: cosh a = r/2a = a, where a is the reduced separation. 
For a longitudinal field, we change the factor from {—jS) to {2(3). The above formula can be 
generalized to account for a pair of dielectric spheres of the different sizes 0. The multiple 
images method was widely adopted p,0,|Tl[ to calculate the total dipole moment as well as 
the interparticle force. 

We should remark that the present multiple images method is an approximation only. In 



fact there is a more complicated images method for a dielectric sphere p]2|JT3|] , which gives 
the exact image dipole moment of a dielectric sphere that placed in front of a point dipole. 
However, the numerical results of the present version of formula agree with the numerical 
solution of an integral equation method |I^ , p!5[ . 



When we apply a sinusoidal electric field, i.e., E(t) = Eosinut, the induced dipole 
moment will vary with time sinusoidally. Due to the nonlinearity of the particles, the 
induced dipole moment will be a superposition of harmonics. In other words, we have: 

Pt = Pw sin tut + sin 3ujt + p^^ sin but + ■ • • (5) 

Due to the inversion symmetry of the dielectric media, only the harmonics of odd order 
survive. The coefficient Pnuj of the nth harmonic is given by: 

Pnu) = — I pTsmnujtdt. (6) 
n Jo 

In what follows, we report results for the transverse field case only. The longitudinal 
field case is similar. In the next section, we will use the series expansion to obtain analytic 



expressions for the harmonics of the induced dipole moment. We will obtain the coefficient 
Pnuj as a power series of the applied field Eq. 



III. SELF-CONSISTENT EVALUATION OF THE AVERAGE FIELD 

According to Eqs.@, the induced dipole moment of the dielectric spheres can be 
determined if we find the average field (E^). The electric field inside the host medium 
can be conveniently calculated by considering the effective nonlinear dielectric constant of a 
two-component composite. For a two-component composite, the effective nonlinear dielectric 
constant is given by ||^: 

where / is the volume fraction of the particles. For a pair of spheres inside a transverse 
field, the effective nonlinear Dielectric constant can be expressed as 

4 = Cm + 1^?^] ■ (8) 

The average field of the host medium can be calculated by using Eq.(|): 

^El) ^ j^EHt)^. (9) 



For a nonlinear characteristic [Eq.(|I])], we solve Eq.(^ self-consistently ||^. The average 
field inside the host medium as well as the dipole moment of the spheres [Eq.(^)] can be 
determined. 

It remains to examine how the higher harmonics of pt depends the nonlinearity. We first 
normalize the material parameters with the linear dielectric coefficient of the host e^: 

c-m 

We then expand (E^) into a Taylor expansion: 

X'JEI) = 4^ + Sl^'-^'^'^ £ {xUeDY , (11) 



where the expansion coefficient Cs is given by: 
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e' =1 
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Hence the expansion coefficients do not depend on the apphed field. Similarly, we expand 
Pi into a Taylor expansion: 
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s=0 



where 
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sinh a 



n=0 



(14) 



sinh(n + l)a 

It is easy to show that: Cg = (s + l)as+i. In the case of a weak nonlinearity, i.e. x^£'^(t) ^ 1 
we can rewrite Eqs.(Pl) and ([l3|), keeping only the lowest orders of x^i?^(t) and Xmi^m)'- 



{El) 



1-f l-f 

The first order term of x'm.E'^it) gives the average field inside a linear dielectric host medium. 
Similarly, the induced dipole moment is given by: 



Pt — aory{t) + ai- 



1-/ ^1-/ 



KiE{t) + K-iE\t) + 



(15) 



It should be remarked that, from Eg. ([131), KiE{t) is the linear dipole moment p^'- 



(16) 



Let us consider a sinusoidal applied electric field E{t) = sin cut. By using the identity 
Asin^ut = 3sinc(j)f: — sinSut, we can expand E^(t) in terms of the first and the third har- 
monics. By comparing Eq.(^ with Eq.(|T3|), we find: 

3 1 
p^ = KiEo + -^K^El and p^^ = --K^E^. 

The above results show that the induced dipole moment must include the higher harmonics 
Furthermore, the results show that p^^ also depends on Eq. This is a nontrivial result 
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as it implies that the first harmonic of the induced dipole moment depends on the strength 
of the nonhnearity. Concomitantly, the higher harmonics should become more significant 
as Eo gets higher. In the case of a higher applied field, we must include even higher order 
terms [i.e., higher powers of E{t)] in the expansion of p^: 

PT = KiE{t) + K-iE\t) + K^E\t) + ■■■. 

Again, by considering the identity 16 sin^ ut = 10 sin ut — 5 sin Sut + sin 5ujt, the harmonics 
are given by: 



-K^E^ + — . 
4 ° 16 

1 o 5 
4 ° 16 



p^ = KiEo + -K:,El + —K^E^ 



= --.K,El - -K,El 



P5u = Y^K^El 

Consequently, in the case of nonlinear ac response, the convergence of the series expan- 
sion is questionable and a self-consistent formalism is needed. The self-consistent equation 
[Eq.(H)] can help to solve for the average field. From Eq.(|Tl]), (|13D and (|^), we find that: 

Pl = Fix'^E',) and xM) = G{x'^E',), (17) 
Po 

where F and G are functions of a single variable, and Pq = ema^PEQ] this is the linear dipole 
moment of an isolated dielectric particle under a time-independent applied field Eq. Similar 
conclusion can be drawn for p^^ and pa^. This demonstrates that we can use x'm^l a 
natural variable in the nonlinear composite problem ||16||. 



IV. NUMERICAL RESULTS 

In this section, we perform numerical calculations to investigate the effects of a nonlinear 
characteristic on the harmonics of the induced dipole moment and the average electric field. 
As shown in section III, the induced dipole moment of the dielectric spheres is affected by 
three factors: the strength of the applied field {Eq), the relative nonlinear dielectric coeffi- 
cient of the host medium (Xm) ^^"^ relative linear dielectric constants of the particles (e^). 
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In order to emphasize the effect of mutual polarization, we use a small reduced separation 
a = 1.1. From Eg. (p!7D , we can use x'm^o variable to plot the numerical results. 

Before showing the numerical results, we consider the simple case of an isolated sphere. 
In this case, the average field inside the host medium can be solved exactly P,p!6|: 

When the applied field is small, the effect of nonlinearity is negligible. The average field is 
dominated by the first part of the above equation. Since the volume fraction / have to be 
small [0, the average field can be well approximated by: 

{El)^E\t). (19) 

For a sinusoidal applied field E{t) = Eosinut, we have: 

Euj ~ Eq, 

which means that the first harmonic varies linearly with Eq. In other words, a linear relation 
between E^^ and Eq is an indication of a weak nonlinearity. 

In Fig.|^, we plot the first and the third harmonics of the average field, with = 2 
(/3 = 1/4) and = 10 (/? = 1/4), respectively. The volume fraction has a small influence 
on the nonlinear effects. The harmonics are plotted against y/x'mEo- The first harmonic 
\/y^^Ei^ varies almost linearly with ^Jx'mEo- Moreover, the magnitude of the third harmonics 
is small compared with the first harmonics. These results show that the average field can 
be approximated by Eg. ([T9|) . 

The numerical results of the average field suggest that the average field varies almost 
linearly with the applied field. However, as suggested by Eq.(|T5|),the dipole moment does 
not vary linearly with the applied field. In Fig.|^, we plot the first and the third harmonics of 
the induced dipole moment. We divide the harmonics by po = ^mO-^PEo, which is the linear 
dipole moment of an isolated dielectric particle under a time-independent applied field Eq. 
In Fig.|^, we find that the first harmonics changes from positive to negative. This is because 
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as the magnitude of the apphed field increase, the field-dependent dipolar factor (3 changes 
from positive to negative. Hence the first harmonics changes its sign. 

Next we examine the third harmonics of the dipole moment. The ratio "Pzu^Ivq increase 
with the applied field. The third harmonics is an evidence of nonlinear effects. The ratio 
increase significantly when the dielectric contrast between the particles and the host is small. 
In other words, a smaller dielectric contrast gives a larger nonlinear response. 



DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. In the present study, we have examined 
the case of linear particles suspending in a nonlinear host. The case of nonlinear particles 



suspending in a linear host has already been studied by similar formalism ||17 . 

So far, we have not considered the frequency dependence of the dielectric constants of the 
host medium and the particles. In a realistic situation, the dielectric constants can decrease 
with the increase of the frequency. For simplicity, we may adopt the Debye relaxation 
expression. 

As we have obtained the expression for the induced dipole moments, we may take a step 
forward to calculate the interparticle force via the energy approach 0. We will find a time- 
dependent force F{f) but only its time average could be measured during an experiment. 
From the energy approach P|, 



{Frit)) 



dr 



We believe that the interparticle force in the etc CclSG should differ significantly from that 
of the dc case of nonlinear ER fiuids because, as we have shown, the nonlinearity enters 
into the composite problem in a nontrivial way. Results of the interparticle force will be 
published elsewhere [jl8| . 



In conclusion, we have considered the effects of a nonlinear characteristic on the ER 
fluid under the influence of a sinusoidal applied fleld. We have calculated the harmonic 
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components of the induced dipole moment as well as the average electric field. We have also 
examined the conditions for obtaining large ac response in ER fluids. 
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FIGURES 

FIG. 1. The first and the third harmonics of the average field, with = 2 (/5 = 1/4) and 
ep = 10 (/3 = 1/4), respectively. The volume fraction has a small influence on the nonlinear effects. 

FIG. 2. The first and the third harmonics of the induced dipole moment. The material param- 
eters are chosen as the same as in Fig.|l|. The nonlinear effects are significant when the dielectric 
contrast is small. 
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